CHAPTER (1 

) : NUMERICAL ANALYSIS 

Curve Fitting 


Method 1 : Least Mean Square 
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Root Mean Square Error ( RMSE ) 
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Some Common Ideas 


(1) y — ax b 

-y 
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(4) y = 


a+fr cos(0 ) 


= Ina + b Inx — Y — A + bX 
1 b 

— — a H — — Y — a + bX 
y x 

X = x 2 

> - — 1 + b cos(0) => - = Y 

y y 


cos(0) = X 


General Form 


OJ^dll ^lt~- jSj 

y = a o 0 0 O) + «101 W + ••• + a n (p n (x) 
y = ^ a 0(x) 

^3 (J (Jl Aju^q 42 <jj£j 42 j^ks-ia <jc- t g W’l j 


■ Z0O 

Z0O01 Z0O02 

Z0O03 _ 

/ a o\ 

/Zy o 0o 

20100 

Z0i 2 

Z0102 

Z0O03 

J 01 - 

/ Zyi0i 

20200 

Z0201 

Z02 2 

Z0203 

r 2 / 

l Zy202 

-20300 

20301 20302 

Z03 J 

W 

\Zy 3 03 




Interpolation : 


1- Lagrange Method 
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Newton Divided method 


P n (x) = y 0 + 8(x — x 0 ) + 8 2 (x — x„)(x — * 1 ) + 5 3 (x — x 0 )(x — x x )(x — x 2 ) 
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Operators 


1- Shift Operator E n /(x) = /(x + nh ) 

2- Forward Operator A/(x) = /(x + h ) — /(x) ; A = E — 1 

3- Backward Operator V/(x) = /(x) — /(x — h) ; V = 1 — E -1 

4- Center Operator 5 /(x) = / (x + — / (x — ; 5 — E^ — E~ 

1 _1 

5- Mean Operator nf(x) = ^ [/ (* + \ ) ~ / (* - fj\ = \ ^ ; 8 = ¥zl± 

A (J.g) = A/.flf r+1 +f.Ag 
A /A = gA/- Ag.f 
\gj 9r+l9r 


operator J gA ^ jSL® 

££ L _1 -1 

A n y 0 = 0 " 2 

^ hVnu L E J) A “3 

(E - 1 ) n y 0 = 0 

1 g 1^*1 (x c_a cJ j *-* j -4 

A n+1 y 0 = 0 jj A n y 0 = 0 ^ j jS -5 






Derivatives by operators form : 




1 1 

A A 2 H — A 3 . 

2 3 

V + ^V 2 + ^V 3 + - 
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by integerating ln(l + x) = x — — + — — 


ln(l — x) = — 
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Newton Forward 
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Newton Backward 


p = 
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Gauss Forward 
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Gauss Backward 
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Stirling Form 
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Approximate for first derivatives 


a) Divided Difference Methods 


Backward divided difference 

Central divided difference 

Forward divided difference 
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Numerical Integration 
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Gaussian Quadrature 
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Three point 
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Romberg 


Therefore 
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Numerical Solution for ODE 




y' = f(x,y ) ; y(x 0 ) = y 0 ;h = — - — 

(one step) 


1- Euler Method 

yt+i =yi + hfi 

Global error 

\E n \=^[{l + hMY-l\ 

2- Taylor Method 
2 nd Order 

h 2 , 

yt+i ~ Vi + hfi + fi 

h 2 h 3 

yi+i = yi + hfi + — fi + — fi 


3 rd Order 




3- 4 th Order RK ( Runge Kutta ) 

Vi+i = y; + g ( K-l + 2 K 2 +2 K 3 + K 4 ) 

Ki = f( x i >Vi ) 

/ h h \ 
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( h h \ 
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4- Modified Euler ( Corrector - Predictor ) 

Predictor Step 

yf+i = Vi + hfi = y£ ( °J 

y[+ [ +1) =yi+\ (fi( x i>yO + f( x i+i 'Vi+i )) 

When you have a 2 nd Order ODE we convert it into a 1 st order system of equations then we 
solve it 
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Numerical Solution of PDE 
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(a ) Fit the curve y = ab 2,x to the following data \ 


Xi 

0.1 

0.2 

0.3 

0.4 

yi ' 

j 0.2 

0.6 

1 

1.4 


and find the root mean squares of errors, 
(b) Use the following table to obtain 
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F(x) 

1.1 
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(i) f(0.01) , f(1.8) and f(0. 55) (use Gauss and Newton’s methods) 

(ii ) f '(0) , f '(0.05), f '(2)and f ”(0.5) h 

(iii ) Find D3,3 (Ricardson extrapolation) whereD nl = a(~— ), 

a(h) is central approximation of f '(1) 

(£> Deduce the recurrence form of Ricardson extrapolation formula 

D nm+1 = ^“D n m - where D u = a(h) = f '00 «■<■ P».i = - fepi) 

11 1 1 »■ ■iiwi |||| " i 1 wi | ' Tmrnmivr* 

With my best wishes ® r: M.Shokry 
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(a )Find P„(x) by Lagrange method for f(x) = cos(rtx) at = 1. 1 , x t 
and X2 = 1. 7 and find the error of calculate at x= 1.15 

( b)From the following table find 
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(i) f(0.01$ f(.75) and f(0.25) (use Gauss forw ard and Newton's methods) 
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(ill) D33 (Ricardson extrapolation) where D14 = f '(0.4). 

(c )Deduce the form of Newton's divided difference low w here 
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With my best wishes 
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